Evaluation of the discovery potential of an underwater 
Mediterranean neutrino telescope taking into account 
the estimated directional resolution and energy of the 

reconstructed tracks 

A. Leisos*, A.G. Tsirigotis, and S.E. Tzamarias, on behalf of the KM3NeT 

consortium 

Physics Laboratory, School of Science and Technology, Hellenic Open University, 
Tsamadou 13-15 and Ag. Andreou, Patras 26222, Greece 



43 

Abstract 

We report on the development of search methods for point-like and extended 
neutrino sources, utilizing the tracking and energy estimation capabilities of 
an underwater, Very Large Volume Neutrino Telescope (VLVnT). We demon- 
strate that the developed techniques offer a significant improvement on the 
telescope's discovery potential. We also present results on the potential of 
the Mediterranean KM3NeT to discover galactic neutrino sources. 

Keywords: Neutrino Telescope, KM3NeT, Galactic Neutrino Sources, 
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1. Introduction 

The geographical position of the Mediterranean is ideal, since the area of 
the sky that is possible to be observed by a deep-sea neutrino telescope will 
include the centre of our galaxy and most of the galactic plane where many 
high-energy gamma-ray sources (potential neutrino emitters) have been dis- 
covered. The optical properties of water and the geographical position could 
elevate a Mediterranean VLVnT, KM3NeT, to be the most sensitive cosmic 
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neutrino telescope. In optimizing the KM3NeT configuration, it is important 
to utilize the maximum of the experimental information in order to maxi- 
mize the discovery potential of the telescope. In this work, we demonstrate 
that by employing in the search strategies information based on the direc- 
tional resolution and the measured energy of each of the detected muons, 
the discovery potential of the telescope is significantly improved. In Section 
2 we describe the detector configuration and we highlight the performance 
of the HOURS [l| software package to simulate the response of the detector 
and to reconstruct the muon-neutrino and muon-antineutrincL] direction from 
the experimental signal. In Section 3 we describe the analysis strategy for 
point-like and extended (disc-type or Gaussian-type) sources. In Section 4 
we present results in estimating the KM3NeT discovery potential for galactic 
sources whilst Section 5 summarizes the results of our study. 

2. Detector Configuration and Signal Reconstruction 

In this study we assumed that the KM3NeT telescope will consist of 308 
vertical structures (Detection Units - DUs), deployed at a Mediterranean 
site (37° north and 18° east) in 3500 m depth. A DU carries photo-sensors 
arranged vertically on 20 horizontal bars (floors), with a distance between 
consecutive floors of 40 m, as described in jij. The DUs fill a cylindrical 
volume of 5.8 km 3 , 1560 m in radius and 760 m in heightH. The positions 
of the DUs on the seabed are almost-uniformly distributed with a typical 
distance between DUs 180 m. Each horizontal bar of a DU, at its end, 
supports two photo-sensor units, hereafter called Digital Optical Modules 
(DOMs). A DOM consists of a 17-inch diameter pressure resistant glass 
sphere housing 31 3-inch photomultiplier (PMT) tubes, their high-voltage 
bases and their interfaces to the data acquisition system with nanosecond 
timing precision. The front-end electronics perform an accurate measurement 
of the PMT-signal arrival time and an estimation of the signal amplitude 
through a Time Over Threshold (TOT) measurement 0]. 

We have simulated the response of this telescope configuration to 15 mil- 



hereafter, we use the term "neutrino" for both muon-neutrinos and muon- 
antineutrinos. 

2 For some special studies in this work we considered the half of the detector, consisting 
of 154 DUs in a cylindrical volume of 2.9 km 3 , with the same DU density as the whole 
telescope. 
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lion neutrinos with directions uniformly distributed on the sky, with energies 
between 100 GeV to 100 PeV, following an energy spectrum proportional to 
E~ 2 . In parallel we simulated, using CORSIKA [4J, the telescope's response 
to muons and muon-bundles produced in 15 million down-going Extensive 
Air Showers (EAS). The simulated EAS have energies between 1 TeV and 
1 EeV, following the experimentally established j^j energy distribution and 
composition of the charged cosmic rays. The direction of the primary cos- 
mic ray on the upper atmosphere is distributed uniformly on the sky up to 
an angle of 3 degrees above the detector horizon. We have used HOURS 
to simulate all the relevant physical processes, the production of signal and 
background, the response of the PMTs as well as the functionality of the 
digitization electronics. HOURS also provides several data analysis tools 
that have been used to reconstruct the direction (estimating also accurately 
the relevant error matrix of the track parameters) of the muons as well as to 
estimate the muon's energy. We have tuned the tracking quality criteria in or- 
der to minimize the probability to miss-reconstruct down-going atmospheric 
muons (or muon bundles) as up-going products of neutrino interactions. The 
method, used in HOURS, to estimate the energy of a reconstructed muon, is 
described in p. For this telescope configuration, the log-energy resolution, 
A (logE), is described by the empirical formula 

A (log E) = 0.32 + 0.082 x arctan (10 - 2.42 x log E) 

where E is the true muon energy in GeV at the closest distance of approach 
to the center of the detector. The log-energy resolution varies from 0.42 at 
1 TeV to 0.23 at 100 TeV, reaching the value 0.2 at 1 PeV. 

3. Analysis Strategy 

As in [7( we choose, for each reconstructed track, a suitable coordinate 
system (hereafter Track Reference Frame - TRF) in which the direction of the 
parent-neutrino with respect to the reconstructed muon direction is described 
by the uncorrelated orthogonal coordinates 9 and <fi. In this system the 
angle if) between the parent-neutrino and the reconstructed muon is very 
well approximated (for small values of if)) as ip 2 = 9 2 + 2 . 

Restricting the analysis to events where the reconstructed muon deviates 
less than -R max from the parent-neutrino direction (i.e. 9 2 + <fi 2 < -R^ax) the 
probability density function (PDF) of 9 and <fi has the following Gaussian 
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form: 





(1) 




where the normalization factor 3ft is well approximated as 



ft= i_ e ^ 




The angle between the parent-neutrino and the reconstructed muon, is 
due to the kinematics of the parent-neutrino interaction as well as to the 
limited detector resolution in estimating the muon track. Consequently, the 
values of the sigma parameters, and E^, in (jTJ, are evaluated as convo- 
lutions of the error matrix elements oq and 0$ of the reconstructed muon 
track due to the detector resolution, with the RMS of the misalignments in 
9 and of the parent-neutrino with the true product-muon momentum, sq 
and respectively (i.e. Eg = + s 2 e and E| = cr^ + s| ). However, due 
to the symmetry of the parent-neutrino direction around the muon track, 
it follows that sq = = s and since the only available experimental in- 
formation relative to the parent-neutrino energy is the estimated energy 
of the reconstructed muon, E m , we used the simulated set of events de- 
scribed in Section 2, to parametrize this misalignment as a function of E m 
(s e (E m ) = (E m ) = s (E rn )). 

When a reconstructed muon track points inside a ring of angular radius 
-Rmax around a point source, Eq. ([1]) expresses the signal-PDF of the track to 
be produced from a neutrino emitted by the source. Restricting the analysis 
only to up-going reconstructed muons and having eliminated the background 
due to miss-reconstructed atmospheric muons, the only background source is 
the atmospheric neutrinos. For small values of -R max EI the PDF describing the 
directional distribution of the background muons around the source direction 
is constant over the solid angle: 



In general, the energy spectra of signal and background neutrinos are 
different, resulting thus in different measured-energy distributions of the 



3 i? maK is chosen large enough to include almost all the muon tracks from the point 
source e.g. of the order of few degrees. 




(2) 
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corresponding muons. We have used samples of simulated events, properly 
reweighted according to the energy-zenith differential distribution of atmo- 
spheric neutrinos [sj, in order to evaluate the distribution of measured ener- 
gies of the background in bins of the cosine of the zenittQ angle, 9 Z . Similarly, 
by reweighting the simulated events, we were able to evaluate the distribution 
of measured energies of muons produced by neutrinos from a source, which 
is located at any declination 5 and emits neutrinos with a certain energy 
spectrum. In order to take advantage of both, the good pointing resolution 
and the energy information provided by the underwater telescope, we define 
the signal PDF, P s (9, 0, E m ; 7), and the background PDF, P b (E m ), as: 

P s (9, 0, E mi 7) = P arig (9, 0) x (E m ; 7, cosfl*) (3) 

P b (E m ) = n (P max ) x P b n (E m ; cos9 z ) (4) 

where P e s n (E m ; 7, cos# 2 ) is the signal PDF of the measured-energy, evaluated 
in 20 bins of cos^ (between -1 and 1) by reweighting the simulated events 
to correspond to a E~~ f neutrino spectrum, and P b n (E m ; cos9 z ) is the back- 
ground PDF of the measured-energy, evaluated in the same bins of cos6* z . 
Furthermore, assuming that a s is the fraction of the signal events in the 
sample of N selected muon tracks, we can write the Likelihood function as: 

l (« s , 7) = n k x ps <?> e ^ 7) ( 5 ) 

i=l,N 

+ (l-a s )xP b (E* m )} 

which can be used to estimate both the signal fraction and the spectral index 
of the source. As a demonstration, we performed 1000 experiments using sim- 
ulated set of events and estimating simultaneously the a s and 7 parameters 
by maximizing the logarithm of Eq. (5) in each set. Each data set consisted: 
a) of 15, on average, background reconstructed muons with directions uni- 
formly distributed inside an angular ring of radius R max = 0.6° around a 
hypothetical point-like neutrino source at 5 = —60° and b) of 15 more re- 
constructed muons produced by neutrinos from the source with an energy 
spectrum proportional to E~ 2 (i.e. 7 = 2). The estimations are presented in 
Fig. [U demonstrating a very good agreement with the true parameter values. 
As discussed in [9], the likelihood estimator, expressed as in (5), can result 



4 The zenith angle refers to to the detector reference frame where Z is parallel to the 
zenith axis and the XY plane is parallel to the horizon. 
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Figure 1: Results of 1000 simulated experiments (see text). In each experiment the number 
of the signal events (N s ) and the spectral index (7) of the source were simultaneously 
estimated, (a) Two-dimensional distribution of the estimations N s and 7. (b and c) 
Distributions of the estimated values of N s and 7 respectively (data points) fitted by 
gaussian functions (solid lines). The mean (n) and sigma (cr) values of the gaussians are 
also shown. 
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in negative values of a Jj in cases that the background is much higher than 
the signal. We have used the expected number of background events, B, in 
the sampled to define the following extended likelihood function: 

u N x e _Al 

£ cx t (a s , 7) = ^ x L (a s , 7) (6) 

where \i = Bj (1 — a s ) is the expected number of total events in the sample, 
a s is the signal fraction and B the mean number of background events. The 
extended likelihood function gives systematically a slightly better esti- 
mation resolution than (5) and the fit procedure converges faster. Following 
the above definitions, \i must be greater or equal to B and consequently a s 
varies in the interval [0, 1]. 

Using ([H]), a resolution parameter A is defined as the following log-likelihood 
ratio: 

A = -2xln /7l (as : 0) (7) 

-^ext \&si 7 'fexp) 

where the numerator corresponds to the hypothesis that all the tracks in the 
selected sample of events are products of atmospheric-neutrino interactions 
and the denominator is the maximum likelihood value under the hypothesis 
that the energy distribution of the reconstructed muons corresponds to a 
known (expected) energy spectrum of the parent-neutrinos (i.e. in the fit we 
estimate only a s ). The discovery potential is defined as the minimum sig- 
nal flux that has 50% probability to be detected during a certain observation 
time (usually one year) with at least 5 sigma significance (p-value=2.86x 10 -7 
in the A distribution) above the background. For any source, the discovery 
potential of a neutrino telescope is evaluated by utilizing Monte Carlo exper- 
iments and comparing the distribution of the A parameter corresponding to 
background-only events with distributions of A corresponding to sets com- 
prising signal and background as described in 0,03. 



In the rest of this Section we demonstrate the performance of this tech- 
nique (hereafter fit-technique) to search for point and extended neutrino 



5 Hereafter, we use the symbol " over a variable to indicate the result of an estimation 
(e.g. x is an estimation of the variable x). 

6 This expected background can be evaluated accurately with simulation studies and 
even more accurately by reshuffling the azimuth of the reconstructed muons in real data 
as in [loj]. 



7 



sources with half of the neutrino telescope (i.e. 154 DUs) described in Sec- 
tion 2. In the case of a hypothetical point source, located at 5 = —60° 
that emits neutrinos with an energy spectrum proportional to E~ 2 and using 
a search ring of -R max = 0.6°, the 5cr discovery potential was found to be 
2.80 x 10~ 9 (E/GeVy 2 GeV _1 s _1 cm~ 2 . We have also ap plie d to the same 
samples of simulated events the so called bin-technique [2, 11], which counts 
tracks pointing inside an angular ring and compares their number with the 
expected number of background muons, without taking into account the 
pointing resolution and the measured energy. The bin-technique, after opti- 
mizing the radius of the search ring, leads to an inferior discovery potential: 
3.80 x 10- 9 (E/GeV)~ 2 GeV-V^m" 2 . 

The gain in sensitivity (observed also in 10] ) is a result of the extra in- 
formation, on the tracking and energy resolution of the detector, which is 
included on a track by track basis. It must be also noticed that the results 
of the bin-method strongly depend on the size of the search ring whilst the 
fit-technique is almost insensitive to such a choice, as long as -R max is greater 
than 2 or 3 times the median of the angle between the parent-neutrino and 
the reconstructed muon direction. However, there are variations of the bin- 
technique where the discovery potential is evaluated by tuning simultaneously 
the reconstruction (e.g. criteria concerning the selection of the "best track 
solution") and the track selection criteria (blindly) in order to achieve opti- 
mum results. Although, this blind optimization leads in some cases to similar 
results as the fit-technique (e.g. the blindly optimized bin-technique finds for 
this source and detector configuration 2.85 x 10~ 9 (E/GeV) 2 GeV _1 s _1 
cm -2 ) such a method is not bias- free when applied to experimental data. 

In order to quantify the contribution of the energy information to the sen- 
sitivity of the fit-technique, we evaluated the discovery potential by employ- 
ing only the tracking information to form the extended likelihood function, 
i.e. P s (9, 0; 7) = Pf ng (6',0) and P b = II(_R max ). In this case the discovery 
potential is found to be 3.20 x lO" 9 (E/GeVy GeV-V 1 cm 2 which when 
compared with the above results indicates that 40% of the gain in sensitivity 
offered by the fit-technique (for E~ 2 neutrino spectrum) is due to the energy 
information. 

Most of the galactic neutrino sources candidates are extended objects. 
In case that the source has a Gaussian shape as several of the H.E.S.S. 
sources (e.g. HESSJ1616-508 or HESSJ1614-518) H, the fit-technique is 
easily modified to account for the source shape by incorporating the source's 
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Gaussian sigma, p, to Eq. (JTJ as: E 2 , = a 2 , + s 2 (E m ) + p 2 and S| = a\ + 
s 2 (i? m ) + p 2 . However, in the case that the source is a uniformly emitting 
disc of radius d, Eq. ([I]) should be integrated over the whole area of this disc. 

Due to the fact that such an integration is very computer-time consuming, 
an independent set of simulated neutrino-induced muon tracks was used to 
parametrize the expected angular profile of the reconstructed muons around 
the center of the disc. The angular signal-PDF, P* ng (9, (f>) is replaced in the 
extended likelihood function (Eqs. [2j |3]and|6]), by ivroffle ( r /^j cos #z)> which 
expresses the PDF of a reconstructed (signal) muon to deviate by an angle 
r from the direction of the center of the source-disc (of radius d), when the 
zenith angle of the reconstructed muon is 6 Z . We used the fit-technique with 
this modification to evaluate the discovery potential of the half detector for 
an extended source at 5 = —60° with radius d = 0.6°, emitting neutrinos 
with an energy spectrum proportional to E~ 2 . Using a search ring with 
-Rmax = 1-2° around the center of the disc, we found the discovery potential 
to be: 4.65 x KT 9 (£/GeV)~ 2 Ge^s^cm" 2 . 

This result is by 67% worse than the discovery potential corresponding to a 
point source at the same declination, reflecting the fact that the source exten- 
sion weakens the impact of the angular resolution to the sensitivity of the de- 
tector for such sources. However, even in the case of this extended source, the 
performance of the bin-technique (5.60 x 10~ 9 (E/GeV)~ 2 GeV _1 s _1 cm~ 2 ) is 
inferior compared to the results of the fit-technique. 



4. Discovery Potential for Galactic Sources 

In this Section we study the potential of the whole KM3NeT detector 
(308 DUs) to discover neutrino emission from galactic objects that have been 
identified as sources of TeV gamma rays. We focus on the bright gamma ray 



source RXJ1713. 7-3946 [13j, assuming that the observed Very High Energy 



(VHE) gamma rays emitted from this object are of hadronic origin 14J and 



that the expected neutrino flux (hereafter reference flux) is described as: 

1.68 x 10- 11 S- 1 - 73 e - V^/" (TeV^cm-V 1 ) (where E is in TeV). We mod- 
eled this source as a radiating disc of angular radius of d = 0.6° centered at 
5 = — 39°46' and we evaluated the discovery potential by the fit-technique, 
reweighting the simulated events according to the reference energy spectrum. 
The 5<7 discovery potential for one year observation was found to be 3.75 
times the reference flux when the energy information is included in the like- 
lihood function (3.9 times the reference flux without including the energy 
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Table 1: Estimation of the discovery potential of KM3NeT and of the required observation 
times for 5a and 3cx discoveries of candidate galactic sources. The discovery potentials are 
estimated for one year of observation and they are expressed as multiples of the reference 
flux of each source. The discovery times are in years and the errors reflect the estimated 
statistical uncertainties. 



Candidate 


Reference Flux 


5a (1 yr) 


5(7 


3<x (1 yr) 


3<T 


Source 


(TeV- 1 cm- 2 s- 1 ) 


Discovery 


Discovery 


Discovery 


Discovery 






Potential 


Time 


Potential 


Time 


RXJ1713.7-3946 


1.68 x io-"£;- 1 - 72 e -V-E/2.i 


3.75 


7.5±1.0 


2.14 


3.0±0.5 


RXJ0852.0-4622 


1.676 x i -a E -i-rs e -^/E/i.i9 


6.30 


33±4 


3.71 


11.0±0.5 


HESSJ1614-518 


0.46 x io-"£;-2-2 e -v , -B/2.i 


16.15 




8.85 




HESSJ1616-508 


0.42 x lO" 11 £;-i.72 e -V £ / 21 


9.54 




5.11 


28±4 



information in the fit). The bin-technique, optimized for best results, offers 



an inferior discovery potential, 4.6 times the reference flux 11]. However, it 



is important to estimate the observation time required to achieve a 5a discov- 
ery (hereafter 5<r-Discovery Time, 5a — DT) in the case that RXJ1713. 7-3946 
emits neutrinos according to the expected flux. Using the fit-technique we 
evaluated the 5a (3a) discovery potentials for several observation times rang- 
ing from 1 to 15 years. The 5a — DT (3a — DT) is the observation time that 
corresponds to a discovery potential equal to the reference flux. When the 
energy information is included the 5a — DT is 7.5 yr, whilst the 5a — DT 
increases to 8 yr in the case that only the angular information is used in 
the fit (in comparison the bin-technique, blindly optimized for best results, 
offers a 5a — DT of 13.1 yr). We have also found that the observation time 
required for a 3a discovery (3<j — DT) with this detector configuration is 3 
yr (with energy information included) and 4 yr (without the energy informa- 
tion). The statistical errors on the above estimates of the discovery time are 
of the order of 1 yr for the 5a — DT and 0.5 yr for the 3a — DT. 

We have applied the fit-technique to evaluate the potential of the de- 
tector to discover other fainter (potential) neutrino sources, which have 
been observed to emit VHE gamma rays, assuming that the observed TeV 
gamma spectrum is of hadronic origin. For example we studied the case of 
RXJ0852.0-4622[l5| (disc with a radius of 1°), HESSJ1614-518 [16] (Gaus- 
sian shape with a = 0.25°) and HESSJ1616-508 [l6| (Gaussian shape with 
a = 0.16°). The prediction of 17J has been used as a reference flux for 
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RXJ0852. 0-4622. For the other two sources we imposed an energy cutoff (the 
same as RXJ1713. 7-3946) to the neutrino spectra by fitting the correspond- 
ing predictions of [17j (which is of the form E~"< due to the limited statistics 
of the observed multi-TeV gamma rays) to the general form N x E~ s e~V ' E / 2A 
with free parameters N and s for each source. We have found the discovery 
potentials of this detector configuration for one year of observation time as 
well as the discovery times in years. The results are summarized in Table 1. 



5. Conclusions 

We have applied the fit-technique in order to estimate the potential of 
a large underwater Mediterranean neutrino telescope to discover point-like, 
(extended) Gaussian-like and (extended) disc-like neutrino sources. We have 
demonstrated that the pointing and energy information offered by this detec- 
tor, when taken into account on a track-by-track basis, enhances significantly 
the discovery potential of the telescope. We have estimated the 5<r discovery 
potential of the telescope for the neutrino emission from the Supernova rem- 
nant RXJ1713. 7-3946 under the assumption that the observed VHE gamma 
spectrum of this source is of hadronic origin. The required flux for 5a dis- 
covery in 1 yr of data taking is: 

3.75 x {1.68 x i -ii£-i-72 e -V^Ai (TeV^cm-V 1 )}. 



We have also estimated that 7.5 yr of observation are required to es- 
tablish a 5<T discovery of this source (in comparison with 13. ly in case the 
bin-technique is used). We have also studied the potential of this telescope 
to discover other candidate neutrino sources in the galactic plane and we 
found that several decades of observation are required. However, it must 
be emphasized that in this study we have imposed an energy cut-off to the 
expected neutrino fluxes from HESSJ1614-518 and HESSJ1616-508 sources 
similar to the RXJ1 713. 7-3946 spectrum. In case of a harder energy spec- 
trum the discovery potential of this KM3NeT configuration to discover these 
sources is better. Furthermore, the discovery potential for galactic sources 
can be improved significantly by optimizing the geometrical parameters of 
the telescope's layout (e.g. the distance between DUs) as well as developing 
more efficient tracking algorithms that take into account the known source 



direction into the signal selection and track reconstruction [18 
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